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I show using Landau theory that quenched dislocations can facilitate the supersolid (SS) to normal 
solid (NS) transition, making it possible for the transition to occur even if it does not in a dislocation- 
free crystal. I make detailed predictions for the dependence of the SS to NS transition temperature 
Tc(L), superfluid density and dislocation spacing L, all of which can be tested against experiments. 
The results should also be applicable to an enormous variety of other systems, including, e.g., 
ferromagnets. 



Recent reports'!^ of supersolidity - a crystal exhibit- 
ing "off-diagonal long-range order" (ODLRO) 2, 3]-in 
solid ^He raise many questions. First, quantum Monte 
Carlo simulations 3 find no supersolid phase. Second, 
the temperature (T) dependence [1,] of the superfluid den- 
sity ps{T) in the supersohd (SS) differs from that in the 
superfluid (SF), contradicting theory Third, no spe- 
cific heat anomaly is seen at the SS to NS transition. 

In this paper, I propose a resolution of these puzzles. 
Since, depending on the material, either local compres- 
sion or local dilation increase the local transition temper- 
ature rc(r)0, and since edge dislocations have regions 
of both types near their coresQ, these defects induce, 
in all materials, regions of elevated Tc, as first noted for 
superconductors 0. ODLRO therefore happens at higher 
temperatures on the tangled network of quenched dislo- 
cations in *He crystals than in the bulk, as in supercon- 
ductors 0, [1] , and can occur even if the clean (disloca- 
tionless) lattice remains normal down to T = 0. 

Specifically, the DGT[3] model with quenched dislo- 
cations implies the following scenario: as temperature T 
decreases below what Fll call the "condensation" temper- 
ature Tcond, which is always > T^'™"^ the transition tem- 
perature of the clean (i.e., dislocationless) lattice, each 
dislocation line in a tangled network of them nucleates 
a cylindrical supersolid "tube" tangent to it. The radius 
of these tubes grows with decreasing temperature. 

We can think of places where dislocations cross, mak- 
ing supersolid tubes overlap, as the "sites" of a random 
lattice. The sections of tube between these "sites" act 
as ferromagnetic "bonds". The typical length of these 
bonds is L, the mean dislocation spacing, which grows 
with annealing; L — > oo for a clean crystal. This random 
lattice does not develop macroscopic supersolidity (or un- 
dergo any phase transition) at Tcond, because the sites 
lack long-range phase coherence near Tcond- However, as 
temperature is lowered further, such coherence inevitably 
develops at T = Tc{L), with Tcond > Tc{L) > Tf"'''^. 
Indeed, if condensation occurs, long-range order always 
develops (i.e., Tc{L) > 0), even if the clean system never 
orders! This ordering at Tc{L) is the SS to NS transition. 

This picture is very similar to Shevchenko'sQ. 

Figure 1 plots the superfluid density Ps{T). When 



j^ciean > g, near Tc{L), 



(1) 



where v « 2/3 is the 3dXY correlation length 
exponent [lot, x = 2(^5^) ~ ^, A and Tq are indepen- 
dent of L, a is a lattice constant, and 



j^clean > (3) 



When temperature T is lowered into the range 

6T{L) <C T - Tf^""" < Tc{L) - rf"'"', (3) 
where ST{L) oc 



'^clean\L^ — 2 



psiT, L) ^ A' T^—-)' " ^ _ r^cleanyl^ (4) 
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FIG. 1: The superfluid density versus temperature for a dis- 
located solid in which the clean system does (top curve) and 
does not (bottom curve) have a transition, ps obeys eqns. 
([1]) and Q in region (I) and (II) respectively, where (II) is 
defined by eqn. ©. Tc(L) for the cases Tf"'''^ > and < 
are, respectively, denoted in this figure by (IS) and TJ^ (L) , 
and given by eqns. Q and ^ . 
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and A' is an L- independent constant. In the L ^ oo 
limit, 6T{L) oc i < Tc(L) - T^='<=°" oc L'^ « L"! , 
ensuring a large window of validity for eqn. Q. Once 
T < Tf"'''^ -5T{L), the tubes overlap, the entire volume 
becomes supersolid, and ps is that of the clean system, 
completely independent of L, and so obeys 



Ps{T) (X {T, 



clean 



(5) 



Note that the high temperature (T > T^'^''" - 5T{L)) 
behavior of Ps{T, L) is strongly sample and annealing de- 
pendent (because L-dependent), but the /oii^temperature 
(T < T^'"^"" — ST{L)) behavior is sample and annealing 
independent, and identical to that of a clean sample. 

Precisely such behavior was recently reported In 

figure 2, PsiT) data from Chan's group [I4I is plotted in 

_ 3 

the form ps * versus T, which Eqn. (j4]) predicts should 
give a straight line section, for T satisfying eqn. ^ . The 
data does indeed show such a straight section, although 
it is fairly short, and the error bars in this region are 
large. More accurate measurements of ps{T) , and of the 
dislocation spacing L (by, e^., ultrasonic velocity and 
attenuation measurements [13|) are clearly needed. Al- 
ternatively, one could deduce the ratio of L's in different 
samples by comparing the coefficients of (T^'^"" — T)~3 
in eqn. (jl]), and using this ratio to test the predicted L 
dependence of ps{T,L) and Tc{L) eqns. H]) and([2]). 
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FIG. 2: Superfluid density versus temperature data from 

. , _ 3 

Chan's group [l^. plotted in the form ps versus T. The 
straight segment of this plot predicted by Eqn. Q is indi- 
cated by the dashed line. 

If the clean system does not order, which I'll refer to 
as T^^'™" < 0, TciL) vanishes as L ^ oo: 



(6) 



where 211 is another L-independent constant, a 
resultref . [9| , could also be tested by measurements 
and/or deductions (as described above) of L. Eqn.([T|) 



still holds near Tc{L), but now with x = 2. Eqns dH) 
and ([5]) never apply, since T = intervenes above T^^''^"". 
The lower curve in Figure 1 plots Ps{T) in this case. 
The experimental situation is currently unclear. 



The Tf 



< scenario is supported by recent 



experiments [14| . showing non-classical rotational inertia 
in unannealed *He crystals, but none after annealing. On 
the other hand, Chan's recent experiments Uj , as dis- 
cussed above and in figure (2), suggest Tf^"-"^ > 0. In 
these experiments of ref. [ll|, single crystals still show 
supersolidity, suggesting that dislocations, rather than 
grain boundaries, are the responsible defects. 

Also suggestive are simulations fis'] which see super- 
solid order near screw dislocations. Although screw dis- 
locations do not, in the DGT model, couple to supersolid 
order, higher order terms allow such coupling([16|). 

The absence of a specific heat anomaly in some ex- 
periments can be explained in this picture. For the case 
rpciean ^ ^.j-^g spccific heat near Tc{L) is given by: 



C oc 



T_ 



Li' 



oc 



Li 



j^clean ^ q 



(7) 



where a = —.0127 is the specific heat exponent of the 
3d XY model (lO|. Clearly, the |r -Tcp" singularity 
vanishes as dislocation density {L oo), and so 
should be seen only in dirty samples, not clean ones. 

The ideas developed here are applicable to, e.g., 
ferromagnets [l7j| , which I'll treat elsewhere 18 1. 

I'll now outline the derivation of these results. My 
Hamiltonian is an isotropic [l^ version of that of 



H = I d^r 



IVV'I 



with 



*(?^ = to + guu{r) 



(8) 



(9) 



Here, to{T) is a decreasing function of temperature T 
satisfying io (T'c = 



< 0, where ft 



is the 



value of io at the transition in the clean system, u and c 
are constants, and ua is the trace of the strain tensor. 

Thermal fluctuations in ua have no effect on the criti- 
cal properties of the superfluid density and specific heat 
at the transition^; I will henceforth ignore them, and 
focus only on strains due to quenched dislocations. 

The clean model will not have a NS —>■ SS transition 
if T^^'^"" < 0. When T^^''^"" > 0, I'll assume (as usual) 



that to{T) = 



where F is a constant, near 



T — X'^c'ean^ p^j. ^ straight edge dislocation running along 
the z-axis with Burgers_yector a along the y-axis, ua — 

, where p and A are the Lame 



4m 



4m 



elastic constantsflQ']. Inserting this into eqn. 



gives 



t(r) ^ to + 



g' cos 9 



(10) 
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where g' ^ ga 

Naively, the system is supersoHd in those regions where 
t(r) < 0. ActuaUy, the mean field transition occurs when 
the minimum energy ■0 (f) first becomes non-zero. The 
temperature at which this occurs is Tcond- 

The Euler-LaGrange equation for eqn. ([5]) is 



, t{r) , u ,n 
c c 



(11) 



As noted in Q ,this equation first has non-trivial (0 ^ 0) 
solutions when io drops below a critical value tcond = 
_2mBo^ where Eq is the quantum mechanical ground 
state energy of a particle of mass m moving in the 2d 
dipole potential V{r) = with P = Varia- 

tional treatments 0, [131, show that Eq — — 7^^|1- where 
0.24 < 7 < 2. So a single dislocation line will, in 
mean field theory, order once io < tcond — Us- 

T — T^ \ : T-j rnclean i 



ing to = r 

7g^^ I rpclean I 
2rc I c I 



|Tcle„„| 

\ rriclean 
^ J- ^ . 



this imphes Tco„d = 
and Tcond > , even if T^^''^"" < 0, 



if ^^|r7 > 1. Hence, condensation onto dislocations can 
happen, even when the clean system does not order. 

However, a one-dimensional system like a single dis- 
location line cannot order. To order, these Id "tubes" 
must cross-link into a three-dimensional network. The 
typical tube length is L, the inter-dislocation distance. 

On length scales ^ the tube radius adt), but L, the 
only important variable is "Goldstone mode"; i.e., the 
phase 6{f) of 0(r) = \^{r) \ e^^^'^K In the tube between 
crosslink sites i and j, 9, on these long length scales, 
depends only on distance s along the tube. This leads to 
a Id Hamiltonian for this tube: 

Hid{{0{s)})^Kid{T) f {dsefds . (12) 

From this, I can obtain an effective Hamiltonian 
Heff{9i,0j) coupling the 9's on sites i and j by inte- 
grating out the 9's along the tube: 



OO „ 



-l3Hia({6(s)}) 



(13) 



where the functional integral D9{s) on the right hand 
side is taken with 9{s) satisfying the boundary conditions 
9{0) = 9i , 9{L) = 9j + 2TTn, where the summation integer 
n in eqn. ([T3|) reflects the 27r periodicity in 9. 

Each of the functional integrals /^^ D9{s) in eqn. (fT3|) 
can most easily be done by rewriting 9{s) as follows: 



9is) 



-9i + 2?™ 
L 



s + 5e[s), (14) 



where the new integration variable 59{s) satisfies the 
boundary conditions 69{Q) = 59{L) — 0. This gives 



DS9is)e-''''^^Io''^''^''^'^^\l5) 



The / DS9 in eqn. p5|) is independent of 9i , 9j and n 
(since the boundary conditions on 69 are) , and so is only 
an overall multiplicative constant in e~^^''ff , which only 
adds an irrelevant constant C to Hp ff{9i, 9i). Hence, 
Hcf f{9i,9j) becomes a "periodic Gaussian" [21| 

Heff{9i, 9j) = Vv{9i — 9j] J) 
= -kBTln[ e-^^'^-'^'-^^'^A+C , (16) 



with the "Villain" coupling 



(17) 



Adding up Hcff{di,9j) for all of the bonds gives a 
model for all of the "sites" (cross links of tubes ): 



Heffim) 



E 

bonds 



vde,--9j-j) 



(18) 



Although these couplings J will be random, due to the 
random bond lengths of the tubes, such "random T^" 
disorder is irrelevant in the RG[2^. and can be ignored. 

This Villain model (|18p orders at a temperature Tc = 
0{J/kB); I will now use this to determine Tc{L). 

Consider first T^^''^"" < 0. In this case, provided 
Tcond > 0, so that Kid{T) ^ 0, we can, for L — > oo, 
estimate Tc by replacing Kid{T) in eqn. (fT7|) with its fi- 
nite, non-zero, T = value Kid{T = 0) = Kq. This gives 
eqn. ® with Tq = j^. Note that taking Kid{T) -> Kq 
in eqn. (|17p is valid since Tc{L oo) 0. 

For the case Tf"^"-'^ > 0, the radu adT) of the tubes 
of supersolid diverge as T ^ j,ciean ^ note 
that the locus on which t(rj_) eqn. (fTO|l is equal to tc 



g' 



which, for to > tc and p > (<) is a 
circle passing through the origin, centered on the negative 
(positive) X-axis of radius 



ac{T) 



9' 



2(io - tc) 



(19) 



Inside this circle, t(f±) < 0, so, naively, this boundary 
(HH) defines the supersolid tube. As T ^ T^'=''=''" from 
above, to tc and so ac{T) diverges: ac{T) oc j_ jdean ■ 
Of course, this argument ignores the V'^ip term in eqn. 
(fTT|) . However, since adT) ^ oo as T ^ j-ciean^ ^ yg^j.-gg 
slowly in space, and we can neglect the term in eqn. 
PT|) and simply balance the other two terms. 

We can include fluctuations in this "local equilibrium" 
approximation simply by replacing the local superfluid 
density ps (f) by its value in a uniform system whose value 
of t equals the local t(r), provided adT) » ^(T) oc 
^rp_-j^cieany^^ whcrc ^{T) and 1/ « § are the correlation 
length and its critical exponent in the the clean system. 
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Since u < 1, ac(T) eqn. ([19]) is indeed >> ^{t) as T 
rpciean fj-Q,^ above. This implies that the local superfluid 
density Ps{r) for T near, but shghtly above, r^='«°", is 



(20) 



where i? is a constant, and I've used the Josephson rela- 
tion ps oc ^^^[23]. This ps acts as the 3 — d "spin- wave 
stiffness" for the phase 0{r); that is, 



H^d — 



1 



local 



(f) 



V6' 



(21) 



with Kiocai{r) — ^ p''°'^°'\r) ■ In the case of a straight 



edge dislocation, taking t{r) from eqn. (fTO|. Ps(7^ by 
eqn. ([20| and 0(r) to vary only with distance s along the 
dislocation line, the Id spin wave stiffness Kid becomes: 



Kid 



<fr^pT\fi_) 



(22) 



Since t(f±^ is constant on circles of fixed radius a, pass- 
ing through the origin, with their centers on the a:-axis, 
and is given by: t{r±) = 2^ — to, I'll change variables of 
integration in eqn. to a. The area of the interval 

[a, a + da] is the difference 2'Kada between the areas of 
the corresponding circles, so I can rewrite eqn. (|22p as 



Kid{T) 



m 



( — — Stn] ada 



A 



(23) 



where Sti) = to — tc, and C{x 
Since Sto oc T - T^c/ea« 

KidiT) oc (T - T^c/ean-)i.-2_ 

earlier expression (|17p for J, 



3.1842: 



for 1/ = 2/3. 
implies that 



2x + l 

eqn. ([231) 

Using this Kid{T) in my 
and then equating the re- 
sult to fcsT, gives eqn. ^ for Tc{L). 

As T drops further, eventually J{T, L) will be > fc^T. 
This is guaranteed to happen, since T can get within 
roughly 5T{L) oc ^ of T^^'^"" before eqn. ([23]) breaks 
down. Since KidiTf"""'^ + (5T(i)) oc L"^-"; J(r^^''="" -|- 
5r(L)) = ^ oc L^-" ^ 00 as L ^ 00, since < 1. In 
this limit, the phase order on the "sites" of the dislocation 
network is nearly perfect, and the standard relationship 
between the macroscopic (as opposed to the local) ps and 
3d spin wave stiffness implies: 



ps{T,L)^J{T,L)/LxO{-j) 

n 



(24) 



which, using (fTT)) for J{T,L), implies eqn. Stan- 
dard results for the model (|16p . ITS)) with the T and 



L-dependent J found above then gives the behaviors of 
Tc{C),ps{T), and the specific heat C(T) quoted earlier. 
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